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ABSTRACT 

The interstellar medium is highly dynamic and turbulent. However, little or no 
attention has been paid in the literature to the implications that this fact has on the 
validity of at least six common assumptions on the Virial Theorem (VT), which are: 
(i) the only role of turbulent motions within a cloud is to provide support against 
collapse, (ii) the surface terms are negligible compared to the volumetric ones, (iii) 
the gravitational term is a binding source for the clouds since it can be approximated 
by the gravitational energy, (iv) the sign of the second-time derivative of the moment 
of inertia determines whether the cloud is contracting (/ < 0) or expanding (/ > 0), 
(v) interstellar clouds are in Virial Equilibrium (VE), and (vi) Larson's (1981) re- 
lations (mean density-size and velocity dispersion-size) are the observational proof 
that clouds arc in VE. Turbulent, supersonic interstellar clouds cannot fulfill these as- 
sumptions, however, because turbulent fragmentation will induce flux of mass, moment 
and energy between the clouds and their environment, and will favor local collapse 
while may disrupt the clouds within a dynamical timescale. It is argued that, al- 
though the observational and numerical evidence suggests that interstellar clouds are 
not in VE, the so-called "Virial Mass" estimations, which actually should be called 
"energy-equipartition mass" estimations, are good order-of magnitude estimations of 
the actual mass of the clouds just because observational surveys will tend to detect 
interstellar clouds appearing to be close to energy equipartition. Similarly, order of 
magnitude estimations of the energy content of the clouds are reasonable. However, 
since clouds are actually out of VE, as suggested by asymmetrical line profiles, they 
should be transient entities. These results are compatible with observationally-based 
estimations for rapid star formation, and call into question the models for the star 
formation efficiency based on clouds being in VE. 

Key words: ISM: general - clouds ~ kinematics and dynamics - turbulence - stars: 
formation 



1 INTRODUCTION 

Interstellar clouds are thought to be turbulent and super- 
sonic. Their Mach numbers range from a 1 (T ~ 7000- 
8000 K, H I clouds) through 10 (T ~ 10 K, molecular clouds 
forming low-mass stars) to 50 (T ~ 10-50 K, molecular 
clouds forming high- mass stars). Since supersonic turbulent 
motions carry mass and produce large-amplitude density 



fluctuations, turbulent fragmentation^ is expected to occur 
in the interstellar medium. 

A useful tool for describing the overall structure of in- 
terstellar clouds is the scalar Virial Theorem (VT). It is 
obtained by dotting the momentum equation by the posi- 
tion vector and integrating over the volume of interest {^^■ 
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^ Turbulent fragmentation is defined as the proeess through 
which a chaotie velocity field produces a clumpy density 
structure i n the gas wit hin a few dynami c al tim e scales 
(see e. g .. Ivon Weizsack eil ll95J: ISasaol Il973l: IScald [Tosl: 
Padoad Il995t iBallesteros-Paredes et al.l Il999at iKlessen 
2000tlHeitsch et al. ll200lHBallesteros-Parede£ll2004d) . 



et alj 
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Being directly derived from the momentum equation, the 
VT always holds for any parcel of fluid. 

Virial Equilibrium (VE) is a restrictive condition of the 
VT, and it is defined by the condition that the parcel of fluid 
under study has a second time derivative of the moment 
of inertia equal to zero. It has been invoked extensively to 
analyze the stability of interstellar clouds. However, since in 
a trans- or super-sonic turbulent interstellar medium, clouds 
should be redistributing their mass as a consequence of their 
own turbulent motions, it already seems difficult to achieve 
VE in a supersonic, turbulent cloud. 

Other simplifications, such as the assumption that the 
cloud is isolated, or that the surface terms in the VT are 
negligible, are frequently made in many (if not in most) as- 
trophysical studies of the VT. Those simplifications have 
been thought to be applicable to molecular clouds and their 
su bstructure for nearly three decades (e.g., the t extbooks 
bvlSpitzejll'975 IShulll99ll: IStahler fc Pailall2005l: iLeaueuxl 
I2OO5I . and references therein), maybe as a consequence of 
the old idea that in the ISM, "all forces are in bala nce and 
the m edium is motionless, with no net acceleration" iSpitzeil 
Il978l . Chap. 11), in which "the observational evidence" 
seemed to be consistent with the expectation tha t inter- 
stella r "clouds tend toward pressure equilibrium" iSpitzeJ 
[l97i). 

The possible inapplicability of these assumptions has 
been mentioned in passing in some previous papers 
(e.g.. iB^testero^Parede^^jJ Il999al : IShadmehri et all 
120021: lBallesteros-Paredej|2004a^ but no attention has been 
paid in general to its implications. Thus, in the present pa- 
per I discuss in detail the applicability of those assumptions 
for molecular clouds and their cores. In S|5|I write explicitly 
the VT for fluids in its Lagrangian and Eulerian forms. In J^l 
I discuss the six more common assumptions of the VT and 
their validity in a turbulent environment. In f0| I explain 
why even though clouds are not in VE, they appear to be 
in energy equipartition, and argue that asymmetries in the 
line profiles are the evidence for clouds out of VE. Finally, 
in 33] I draw the main conclusions. 



2 THE VIRIAL THEOREM 

The Virial Theorem can be derived from the momentum 
equation, by dotting it by the position vector x and inte- 
grating it over the volume of interest. Althoughit is usually 
written i n its Lagrangian form, i.e., by following the mass 
(see, e.g.. lSpitzeill978HShull99ll:lHartmannll99^ . it can be 
also obtained in its Eulerian form, i.e., by fixing the volume 
in space (see, e.g.. .Parker.. 1979. : .McKee and ZweibcL.1992) . 
obtaining 



1 iflE 

2 dt^ 



Id^lL 



2 Ey, 



-2( Tkin + Tint ] + M + rM ~W {!) 



2 dt^ 



= 2 Bkin + Eint - 2ri 



M + TM - W 



(2) 



where I ~ Jy pr^dV is the moment of inertia of the cloud 
(subindexes E and L in eqs. and [5| stand for Eule- 



rian and Lagrangian, respectively), -Ekin = 1/2 Jypu^dV 
and Tkin = —1/2 J^XipuiUjhjdS are the kinetic energy of 
the cloud and the kinetic stresses evaluated at the sur- 
face of the cloud, respectively, i?int = 3/2 Jy PdV is the 
internal energy. Tint = —1/2 J^ XiPhidS is the pressure 
surface term, M = 1/8tt Jy B^dV is the magnetic energy, 
TM ~ 1/4tv J XiBiBjfijdS is the magnetic stress at the sur- 
face of the cloud, W — Jy Xip d4>/dxi dV is the gravita- 
tional term, with being the gravitational potential, and 
^ ~ §s P^i'''^f''i'^^ is the fiux of moment of inertia through 
the surface of the cloud. In the previous equations, p, Ui, Bi, 
P, and n are the density, the i*^ component of the velocity 
u, the i*'' component of the magnetic field B, the pressure, 
and a unitary vector perpendicular to the surface S that 
surrounds the volume V , over which the integrals are per- 
formed, respectively. In the notation above it is used the 
Einstein convention, where repeated indexes are summed. 



3 THE COMMON ASSUMPTIONS 

Various assumptions for the terms involved in eqs. Q and 
have been made in the literature. Some of them, indeed, 
have been converted into myths, since their applicability to 
interstellar clouds not only is not demonstrated, but it is not 
even questioned, either in textbooks, or research papers. 



3.1 First assumption: The kinetic energy is 
generally a term of support 

It is almost universally considered in the literature that the 
turbulent (or kinetic) energy, Skin = 1/2 / pu^ dV, provides 
support to clouds against collapse. While this is true for a 
system of particles, and partially valid if the kinetic energy 
is in the form of large-scale expansion and/or rotation, it is 
by no means certain that all the kinetic energy available will 
help against collapse in a system where turbulent fragmen- 
tation can occur, as discussed below. 

This idea has its origins in IChandrasekhail il95 J) , who 
proposed that in the analysis of the gravitational instabil- 
ity the turbulent velocity field should be included. In his 
description, an effective sound speed is introduced, given by 

Ccff ~ ~\' ^'^rms (3) 

where Cs is the sound speed, and Umis is th e velocity dis- 
persion of the turbulent motions (see, e.g., iKlessen et alJ 
2000; Mac Lo w and Klossen .2004 . for a review). This de- 
scription is valid only if (a) turbulent motions are con- 
fined to scales much smaller th an the size of the cloud 
iBallesteros-Paredes et alJll999a)) . and (b) such motions do 
not produce new, smaller-scale Jeans-unstable density en- 
hancements. The first hypothesis disregar ds one of the main 
features of tu rbulent flows in general ('e.g.. lKolmogoro v|l94l|: 
lLes ieuill990l) . and of interstellar clouds in particular jLareonl 
198]|), namely, that the largest velocities occur at the largest 
scal es. An attempt to in clude this fact has been proposed by 
iBonazzola et alJ il987r) . who suggested including the value 
of the rms velocity dispersion at each scale I oc 1/k, i.e.. 
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(4) 



where k the wavenumber corresponding to the scale I, and 
{u{k)) is given by the energy spectrum E{k) = Ck~^ as 



{n(fc)) = 



E(k)dk 



C 



:k'- 



(5) 



where C is a constant and 5 is the spectral index. The second 
condition, i.e., that turbulent motions do not produce Jeans- 
unstable density enhancements, has the underlying compli- 
cation that motions at scales larger than / ~ 1/fc will be very 
anisotropic with respect to structures of size I. Those modes 
will produce shear (through vortical modes) or compressions 
(through compressible modes) to the structures^. Compres- 
sions in particular red uce the local Jeans mass (Sasao 1973 
iHunter fc Flec3ll982ri . and can induce local collapse. Thus, 
a fraction of the turbulent kinetic energy is involved in pro- 
moting collapse, rather than opposing it. 

By decomposing the velocity field in its solenoidal and 
compressible components, the kinetic energy modes that 
provide support to the clouds are those having divergence 
larger or equal to zero, 



V • ?i ^ 0. 



(6) 



This includes the solenoidal modes ( V • u = 0) , and the ex- 
pansional component of the compressible modes (V • u > 0). 
In other words, the precise result of collapse or support must 
then reflect the balance between all the agents that favor 
collapse against those agents that provide support. In the 
first group, not only the gravitational energy should be in- 
cluded, but also the kinetic energy involved in the compress- 
ible modes (-E'kin,V ti<o), versus the kinetic energy involved 
in the expansional and rotational modes (-Ekin.Vu^o)- 



3.2 Second assumption: The surface terms are 
negligible 

It is frequently found in the literature that the surface terms 
are neglected a ltogether, especially in ob s ervational work 
jLarsonI Il98ll: iMvers fc GoodmanI Il988t iFuUer fc Mw3 
[^9^? !mainlv because there is not a direct way of measuring 



them observationally, althou gh it is also a common practice 
in th eoretical studies (e.g., ICha ndrasekhar fc Fermj|^^53 ; 

197C; 



^aricor 1969 ) and textbooks (e.g ., Suitzcr 1978; Parke 

Ishulli99ll: IStahler fc Fallal 120051) . This assumption is based 
on the idea that self-gravitating clouds may be considered 
isolated because then their internal energies dominate the 
dynamics. The most notable exception is the thermal pres- 



sure confinement" 


e.e..lMcCroa"l957':'Keto fc Mverdll986l: 


iBertoldi fc McKed 


19921: lYonekura et al.. 


.1997|). AhhouEh 



some works have considered, by analogy, the possibility of 
turbulent pressure confinement by means of the term Tkin 

^ It should be remembered that compressible and v ortical 
modes are coupled, and they exchange ener g y (e.g., ISasaol 



tney 

| 1973t IVazQuez-Semadeni et alj Il996t iKornreich fc Scald " I2OOO 
lElmegreen fc Scaldl20o3) ! 



re.g.. lMcKee and Zweibeil992h . such confinement of a cloud 
is difficult to achieve because the large-scale turbulent mo- 
tions are anisotropic an d will in genera l distort or even dis- 
rupt the cloud (B allesteros-Paredes et al...l999ai . 

Although we cannot measure the surface terms from 
observations, the possibility that they are as important as 
their corresponding volumetric terms suggest to investigate 
into numerical simulati ons of the interstellar medium. In 
fact, Ballcstor os-Paredes fc Vazguez-Semadenil |[l997'l found 
that for an ensemble of clouds in two-dimensional simula- 
tions of the interstellar medium at a kiloparsec scale, the 
surface terms have magnitudes as large as those of the vol- 
umetric ones (for thre e -dimensional simulation s , see also 
I Shadmehri et al.l l2002l : iTiUev fc Pudritd 120041: lOib et alJ 
I2OO6I) . This result suggests that, on one hand, either sur- 
face and volumetric terms are of comparable importance in 
shaping and supporting the clouds. On the other hand, it 
suggests that clouds must be interchanging mass, momen- 
tum and energy with the surrounding medium. In such an 
environment, the meaning of thermal or ram pressure con- 
finement is not clear, since motions at all scales must morph 
and deform the cloud. 



3.3 Third assumption: The gravitational term is 
the gravitational energy 

The gravitational term entering the VT is written as 



W - 



P 1^ dV. 

OXi 



(7) 



Splitting-up the gravitational potential as the contribution 
from the cloud itself (<^cioud), plus the contribution from the 
outside {(pcxt), 

— 0cloud -I- 0oxt , (8) 

the gravitational term can be written as: 



W = E,- 



p — — dV 

OXi 



(9) 



where Eg = —1/2 p(f> dV is the gravitational energy 
of the cloud alone, since the volume of integration of 
^cioud and the volume V of the integral coincide (e.g., 
[Chandrasckhar fc Fermi Il953f) . The second term in the 
right-hand side is usually either implicitly or explicitly as- 
sumed to be negligible compared to the first term, and the 
gravitational term W is then assumed to equal the gravi- 
tational energy Eg. This is valid only if the distribution of 
mass is spheroidal, or if the medium outside the cloud is ten- 
uous such that its contribution to the potential is negligible. 
However, clouds are more similar to irregular fractals with 
arbit rary shapes (frequently filamentary) than to spheroids 
(e.g.. iFalgarone et alJll991ll . and the contribution from the 
external gravitational field may not be negligible, giving rise 
to tidal torques. Although up to now there is no observa- 
tional estimation of the contribution of the external mass to 
the gravitational term for any interstellar cloud, mass esti- 
mates for for H I "envelopes" around molecular clouds are 
of the same order of magnitu de than the mass of the molec- 
ular clouds themselves (e.g., IWilliams fc Maddalenal|l996l : 
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iMoriartv-Schieven et al.lll99^ . Thus, it is not diHicult to 
realize that the contribution of the second term in the right- 
hand side of eq. Q to the gravitational term W can be of 
the same order of magnitude as the gravitational energy Eg. 

Similar arguments can be made for the interiors of 
molecular clouds: even though we can approxim ate their 
shapes as (triaxial) spheroids jjiiina et alJll999l) . embed- 
ded molecular cloud cores are subject not only to their own 
self-gravity, but also to the tidal forces from their parental 
molecular cloud. Thus, it is not clear that the tidal forces 
represented in the second term of eq. @ will be negligible, 
and the assumption that the gravitational term equals the 
gravitational energy of the cloud seems unjustified. 

What is the meaning and the effect of the second term 
of eq. @ on the energy budget of the clouds? It can be seen 
that it involves the gradient of the external potential. For 
non-symmetrical distributions of mass, this term is out of 
balance even if the distribution of mass inside the cloud is 
symmetric. Thus, this term represents the tidal forces over 
the mass contained in the volume V, and it can be split up 
into three terms: 



Xi p d(t>cxt/dxi dV = 9(00x1 Xi p)/dxi dV 



3 / ^cxt p dV ~ ^cxt Xi dp/dxi dV . (10) 



Using the Gauss theorem, the first term on the right- 
hand side of eq. 1101 can be interpreted as the gravitational 
pressure evaluated at the boundary of the cloud. 



9(0cxt P Xt) 

dxi 



dV ^ (p ((^cxt p Xi) lii dS. 



(11) 



3.4 Fourth assumption: the sign of J defines 

whether the cloud is collapsing or expanding 

It is frequently argued in the literature that if / > 0, the 
cloud must be expanding, while / < implies that the cloud 
is contracting. Equilibrium is assumed to occur when 7 = 0. 

This idea rests on the fact that the gravitational energy 
Eg has a negative sign, and it is a confining agent, while the 
internal energies (thermal, kinetic, or magnetic) are positive, 
and they are assumed to act as supporting agents against 
collapse. By neglecting the other terms, it can be assumed 
that if the gravitational energy is larger than the sum of the 
internal energies, the sign of the right-hand side of eq. ||5J is 
negative. Physically, if gravity wins, the cloud collapses. 

Although energetically this is true, it is not hard to 
find an example in which an expanding cloud has a negative 
second time derivative of the moment of inertia. Assume a 
sphere with constant density and fixed mass M. Its moment 
of inertia is 

/ = \-kMR'. (14) 



If its size varies with time, for instance, as a power law, 



R — R{p) — ^0 



to 



its second time derivative is given by 

2 / _ \ 2-1-2 / 

7(27-1 



6M / Ro 



to 



(15) 



(16) 



which is negative if0<7<l/2, even though it is expand- 
ing. In general, I has been treated in the literature as if it 
were I. 



The second term, by similarity with the gravitational 
energy, can be interpreted as three times the work done to 
assemble the density distribution of the cloud against the 
external mass. 



p dV. 



Finally, the last term in eq (|7J, 

/>cxt Xi dV 

axi 



(12) 



(13) 



involves the gradient of the density field inside the cloud. 
Although there is no clear interpretation of this term, it is 
worthwhile noting that its contribution is null for a homo- 
geneous distribution of mass inside the volume V of integra- 
tion. 

From the numerical point of view, several stud 



les l|B allester os-Paredes fc Vazauez-Semadeni 
Ballesteros Pareded |l999l: IShadmehri et al. 



1997 



2002 : 



Tillev fc Pudritj2004l:lDib et alJi200dl have found that the 
gravitational term can be negative or positive. If negative, 
it will be a confining agent. If positive, its overall action 
will be to contribute to the disruption of the cloud/core. 



3.5 Fifth assumption: interstellar clouds are in 
Virial Equilibrium 

The definition of Virial Equilibrium is that the left-hand 
side of the Lagrangian Virial Theorem (eq. 0) equals zero: 



/l =0 



(17) 



(see, e.g.. ISpitzerlll97^ . Although there are some ob serva- 
tiona l papers show i ng molecular cloud s out of VE (e.g.. lCar'il 
Il987t lLorenlll98sl: iHever et alj|200ll) . it is frequently en- 
countered in the literatur e the VE assumption for MCs and 
their substructure (e.g., 'Mvers & Goodman 1988; McKed 
■ 1999; Krumhol z fc McKee ,2005.: .Ward-Thompson et all 
Eood iTan et eoT l2006l) . As mentioned in the introduc- 
tion, this statement is based on the (old) idea that 
all the forces (in the ISM) should be in balance, and 
the medium should have no net acceleration (see e.g., 

iSpitzcr. 1978,. , Cha p. 11) . Howe ver, b o th observational 

_li'jenkins etall Il983l: iBgwvcr ot aP Il995l: Ijenkins fc Trlppl 
'200lHJenkinsl 120021: iRedfield fc Liriskvll2004). and numer- 



ical (IW zml^z^Semad^i^^^JT20 03l: iMac Low et alJ I2OO.- 
iGazol . Wzau^ ^^emad^n^^ ^n7'200^ studies have found 
that the ISM is not in strict pressure balance, but exhibits 
strong pressure fluctuations, and in fact the turbulent ram 
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pressure is siRnificant ly larger than the thermal one (e.g., 
iBoulares fc CoJl99d) . 

From an observational point of view, it cannot be 
demonstrated that clouds are in VE because neither the 
detailed three-dimensional density structure of molecular 
clouds, nor the time derivatives for interstellar clouds can be 
measured observationally. From a theoretical point of view, 
the VE assumption has strong implications: either super- 
sonically turbulent clouds are not redistributing their mass 
inside, or the way in which the time-derivative of the mo- 
ment of inertia (J) vary (which can be interpreted as the 
time- variation of their mass redistribution) is constant. Both 
statements seem implausible in a hig hly dynamical, non- 
linear ISM. iMcKee and ZweibeJ lll993') . for instance, have 
recognized the difficulty of achieving VE in a turbulent 
medium, because turbulent motions carry fluid elements, re- 
distributing their mass. 

From the numerical point of view, fortunately, it is 
possible to calculate all the terms of the VT for clouds, 
since we know all the variables involved in the numer- 
ical simulations. In order to tost the Virial theorem, 
i Ballc stcros-Parcdcs & Va zquez- S em_adc nii_ see also 

iBallesteros-Paredes^etTarj Il999al : iBallesteros ParedesI 119991) 
calculated all the integrals in the EVT for clouds in numer- 
ical sim ulations of th e inter stellar medium at a kiloparsec 
scale by iPassot et al.l (Il995) . They found that the second 
time derivative of the moment of inertia never goes to zero, 
sug gesting that clo uds must be transient entities. 

McKed (1 19991) has recognized the difficulty to achiev- 



ing VE for actual MCs. He proposed two possibilities 
for assuming clouds in VE: the flrst one is to average 
in time the second time derivative of the moment of 
inertia. He suggests that {/) = if the averaging time 
considered is much larger than the dynamical timescale of 
the cloud, i.e., if tavg tdyn. The second one is that for 
an ensemble of clouds, some of them may have positive 
values of I, and some others will have negative values, 
so that VE holds for the ensemble. Regarding the first 
assumption, there is a hidden assumption behind it: the 
cloud maintains its identity during several dynamical 
timescales"^ . However, from the numerical point of view, 
simulations of the ISM suggest that the clouds are continu- 
ally morphing and exchanging mass, energy and momentum 
with their surroundings, so that over t ~ tdyn they have 
changed significantly. On the other hand, observational 
evidence llBallesteros-Par edes et al]|l999l: |Elmegreen||2000|; 
[Hartmam^^dl |20^ [BaUestero^Pare^^&^artoa^ 
I2OO6I) also suggests that the lifetimes of the clouds are not 
signiflcantly larger than their own dynamical times (i.e., 
they are transie nt). Concerning the second assumption, 
the a nalysis bv [Ballesteros-Paredes fc Vazauez-Semade'nil 
il997l see also iBedl^s^ro^PETOdeT^ggf) shows that the 
moment of inertia spans up to 7 orders of magnitude (in 



In fact, it has to be oscillating around a mean shape without 
a strong redistribution of mass, in order to achieve the condition 
(I) = 0. A similar assumption is made by iMcKee and Zweibell 



absolute value) between the largest clouds and the smallest. 
Thus, it is not clear that an average for such large scatter 
will be representative of the actual dynamics of the clouds. 
In other words, even if I did average out to zero for the 
ensemble, this does not alter the fact that, individually, 
clouds are not in VE. 



3.6 Sixth assumption: Larson's relationships are 
the observational demonstration of clouds 
being in Virial Equilibrium 

The mean density -size and veloc ity dispersion-size relations 
first discussed by iLarsorJ il98ll) are thought to be an ob- 
servational demonstration of Virial Equilibrium. However, 
as iMvers fc Goodman ( 1988) recognize, they are actually 
only compatible with equipartition. For instance, assuming 
equipartition between kinetic (turbulent) and gravitational 
energy. 



Sv 



GM 2 

— « {p)R 



(18) 



Thus, if there is a mean density-size power-law relationship 
p oc 7?", then, there should be a velocity dispersion-size 
relationship of the form 



Sv^ oc Ri^, 
with 

/?=(a + 2)/2. 



(19) 



(20) 



Several caveats must be mentioned at this point. First of 
all, this derivation does not mean VE, but energy equipar- 



tition, since the only assumption was Eg 



(see 



eq. 1181 ). A similar result can be found if the equiparti- 
tion assumed is valid between the gravitational and mag- 
netic energy. In this case, 5v is proportio nal to the Alfven 
speed fsee . e.g.. iMvers fc Goodman! Il988l). S econd , as dis- 
cussed by IVazouez^emaden^^GazoJ lll995l) , in the par- 
ticular case of p oc -R"'^, dv oc R^^^. However, the pair 
a — —l,l3 — 1/2 is not unique. Any pair of values 
satisfying the equation 1201 will be consistent with en- 
ergy equipartition (again, not VE). Finally, it is convenient 
to recall that the validity of Larson's relations has been 
called int o question, especially the mean density-size rela- 
tionship iK cgcl '198S^; 'Scalo"l99d'; 'Vazgue z-Semadeni et alJ 
I1997I: IBallesteros- Paredes fc Mac Lowii200j) . It seems that 
this relationship is more a consequence of the observational 
process, in which the dynamical range of the observations is 
limited below by the minimum sensitivity of the telescopes, 
and above by saturation of the detectors, optically thick ef- 
fects, and depletion. 



4 DISCUSSION: VIRIAL MASS, VIRIAL 
EQUILIBRIUM AND THE STAR 
FORMATION EFFICIENCY 

As it is discussed, there are at least six common assumptions 
related to the Virial Theorem which seem to be unjustifled 
for a turbulent ISM. Some questions, such as why Virial 
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Masses are good order-of-magnitude estimations of the ac- 
tual mass, or why in principle, the sub- or super-criticality of 
a molecular cloud core is a good estimation of the dynamical 
state of such a core, are still valid to ask. 

The answer is that those estimations are based on en- 
ergy equipartition, not on VE (i.e., iSkin ^ Eg ^ Em 
does not mean that II is negligible when it is compared 
to the left-hand side of eq. Molecular clouds seem to be 
in approximated equipartition b etween self-gravity, kinetic , 
and magnetic energy (see, e.g.. iMvers fc GoodmanI Il988ll 
although the super-Alfvenic nature of molecular clouds is 
still a matter of debate, (e.g., ^P adoan & Nordlund 1999, 
see alsoJWar^^hon^soi^^Ll200Q for an observational re- 
view; iBallestero^Pare^^^ayli^ for a theoretical one). 

The question thus, is why MCs seem to be in energy 
equipartition? The answer is related to what we identify 
as a cloud, and to observational limitations: in the first 
case, clouds with a substantial excess of internal energy will 
rapidly expand and merge with the more diffuse medium. 
This is the case of an H II region or SN explosion expanding 
within its parental molecular cloud. The amount of energy 
provided by those events is enough to disrupt their pa rental 
environment within a few Myr fe.g.. lFranco et al.ll994ll . The 
excess of internal energy is mostly in the form of UV pho- 
tons, which rapidly ionize the molecular gas, reducing even 
more the life timescale of what it is identified as the parental 
molecular cloud (see, e.g.. iMellema et alj|2005l) . In the sec- 
ond case, if the cloud has an excess of gravitational energy, 
its free-fall time velocity is at most a factor of y/2 the ve- 
locity dispersion needed for equipartition. This means that 
the difference between a free-fall collapsing cloud and one in 
equilibrium is just a factor of two in the energy, and a fac- 
tor of ^/2 in velocity. Both systems will be, in principle, in 
order-of-magnitude energy equipartition. Note, furthermore, 
that in the free-fall collapsing case all the kinetic energy will 
be observationally identified as energy for support, where it 
is not actually providing any support. This re-enforces the 
need to distinguish between the compressible, expansive and 
vortical modes of the kinetic energy (see H3.H . 

As discussed above, not all the terms entering 
the VT can be measured observationally. However, all 
of them can be measured from numerical simulati ons. 
The w ork by rBallcstcros-Parcdcs & Vazg ucz- Scm adenil 
^1995^: jBallcstcros-Parcdcs & Vazgucz-Scma denil lll997f) : 
IShadmehri et al.. (.2002.1 ; Tillev fc Pudritz. (.20o3H or turbu- 
lent realizations of MCs shows that all the terms entering 
the VT are of similar importance. Note that the same prob- 
lems related to the identification of an expanding cloud are 
applicable to the numerical work. In this case, what we iden- 
tify as a cloud are the density enhancements. By definition, 
vacuumed regions (due either to modeled stellar activity, or 
by the turbulence itself) are no longer considered clouds. As 
for the collapsing case, strongly self-gravitating clouds will 
develop a velocity field that is a factor of y/2 the value of en- 
ergy equipartition, if only gravity is acting, and of the order 
of magnitude if the turbulence is forced. Thus, order of mag- 
nitude equipartition between magnetic, kinetic, and gravita- 
tional energy is valid to assume to predict whether the clouds 
should be supercritical or subcritical jBertoldi fc McKed 



119921: In akan3ll99^ . or to understand why subcriti cal cores 
are u nlikely to survive long-lifetimes within MCs llNakand 
llQSM) . However, the order-of-magnitude coincidence of the 
involved energies does not mean that clouds are in equilib- 
rium at all. It should be stressed that the difference is not 
just semantic (Virial mass vs energy-equipartition mass). 
The difference is conceptual: turbulent clouds cannot be in 
equilibrium because turbulent fragmentation takes them out 
of equilibrium, and thus they do not last long. In order to 
achieve equilibrium, (even if it is a time-ave raged equilib- 
rium, (i) = 0, as suggested bv lMcKe3 (|W9^), it is neces- 
sary to allow the cloud to live for, at least, several crossing 
times. To give an example, the Taurus Molecular Cloud has 
a size of i ~ 20 pc, and a velocity dispersion of Sv ^ 2 km 
sec~^. Its lateral crossing time is Tcross = l/Sv ~ 10 Myr. In 
order to achieve equilibrium, Taurus has to live 20-30 Myr 
without a serious distortion or modification. Such a condi- 
tion seems difficult to achieve if one realizes that Taurus has 
a Mach number of the order of 10, and the associated H I 
gas has Mach numbers of 20-30 in respect to the internal 
CO gas (sec, e.g., the velocity-position diagrams of H I and 
CO in Ball cstcros-Parcdes et al. 1999). 

The fact that turbulent clouds cannot be in VE and 
cann ot last seve r al cro ssing times is in clear contradiction 
with iTan et al] (l2006l) . who made theoretical arguments 
that favor star formation occurring during several dynami- 
cal timescales. However, in their theoretical derivation they 
applied the star for mation rate per free-f a ll tim e expression 
given by eq. (30) of iKrumholz fc McKe3 ll2005ll to a clum^ 
in hydrostatic equilibrium. As iKrumholzfcMc Kcc (200. 
point out, their eq. (30) is valid only for large Mach num- 
bers'* (20-40), while by definition, a clump in hydrostatic 
equilibrium is subsonic. Thus, the theoretical arguments in 
favor of slow star formation are wrong 

Another point to be stressed concerns line profiles of 
H I, CO and its isotopes, and even higher density tracers 
(CS, NH3, etc.). If interstellar clouds and/or their cores are 
not in Virial Equilibrium, they must be distorted and dis- 
rupted within a dynamical timescale. How will their line- 
profiles look, and how will they look if they were in equilib- 
rium? Certainly, if turbulence were microscopic (necessary 
to achieve equilibrium), line profiles should be symmetric. 
The fact that H I and CO clouds (e.g., iHartmann fc Burton! 
Il997t Id ame et al."'200l') as well as molecular cloud cores 
(e.g., Falgarone et al. 1998) exhibit non-symmetric line pro- 
files is probably the best and the only evidence that the term 
II is different from zero^. Even a small degree of asymme- 



•* Other problems involve d in the derivation of eq (30) by 
iKrumholz fc McKej J200^ will be discussed in a further con- 
tribution. 

^ Note that in an Eulerian frame of reference, an idealized dust- 
lane in a spiral wave may appear to be long-lived in terms of 
its overall structure even though particular gas molecules do not 
stay inside it very long. For such system the line profile will be 
asymmetric, while its Eulerian second-time derivative of the mo- 
ment of inertia can be zero {Ie = 0). This system, however, is not 
in VE, as the definition of VE, eq. 1171 involves the Lagrangian 
frame of reference (Note that the difference between the Eulerian 
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try in the line profiles suggests that large-scale motions are 
present in the observed system, which is a natural conse- 
quence of turbulence being a multiscale phenomenon. This 
does not mean that such a system is co llapsing or expand- 
ing as a whole, as formerly suggested bv lGoldreich fc KwanI 
l|l974) . Instead, it means that turbulent large-scale motions 
are present in the system, which should be evolving within 
a dynamical timescale. 

In this context, it should be stressed that clouds pre- 
senting large-scale motions and evolving to form stars 
rapidly not necessarily will have a hi gh star formation ef- 
ficienc y, as is the common belief since IZuckerman fc Evari^ 
il974h . Although former models of quasi-static evolution of 
molecular clouds were proposed to reduce the star forma- 
tion efficiency, the turbulent models producing local col- 
lapse rapidly, have small efficiency because gravo-turbulent 
fragmentation involves only a sm all fraction of the meiss of 
the system in collapsing regions (jVazauez-Semadeni et all 
l2003l . [200^. In numerical models, when energy feedback 
fro m stars is included, the clou d is blowed out rapidly (see, 
e.g.. lBallesteros-Paredes^l2004a^ . However, it should be rec- 
ognized that a detailed quantification of the star formation 
efficiency in turbulent simulations with open boundary con- 
ditions and stellar feedback is needed. 



5 CONCLUSIONS 

The present contribution has discussed the applicability of 
the six more common assumptions on the Virial Theorem. 
Specifically, 

(i) It was shown that a decomposition of the velocity field 
into its vortical and compressible modes is necessary, since 
only modes satisfying the condition V • it > provide sup- 
port, while modes satisfying V ■ u < foment collapse. 

(ii) It was argued that for a supersonic, turbulent ISM, 
surface terms should not be neglected. 

(iii) It was shown that the gravitational term can be de- 
composed into a contribution from the cloud itself, and a 
contribution from the outside. The first part is the well- 
known gravitational energy. The second part is the sum of 
three terms: The gravitational pressure evaluated at the sur- 
face of the cloud plus three times the work done against the 
external mass to assemble the density distribution of the 
cloud, plus a term that depends on the gradient of the den- 
sity distribution of the cloud. These represent the tidal forces 
due to an external potential, as could be the case of a dense 
core within a giant molecular cloud, or a giant molecular 
cloud close to a spiral arm. It is argued that this contribu- 
tion can be as important as the self-gravitational energy. 

(iv) Using a simple counter-example, it was shown that 
the sign of the second-time derivative of the moment of in- 
ertia does not determine whether the cloud is contracting 
or expanding. An expanding cloud may very well satisfy the 
condition / < 0, and a contracting one may satisfy J > 0, 

an d Lagrangian VT is Jf, — / p; = — 4r]jin -|- d^/dt, as pointed out 
bv lMcKee and Zweibeilll992l) . 



contrary to the common belief. In other words, I has been 
treated in the literature as if it were I. 

(v) It was argued that interstellar clouds are not likely to 
satisfy the Virial Equilibrium (VE) condition I — 0. 

(vi) It was shown that Larson's (1981) relations are not 
observational proof for clouds being in VE. 

(vii) Clouds seem to be in energy equipartition because 
of either observational limitations, as well as because of the 
intrinsic definition of a cloud. 

Turbulent fragmentation plays a crucial role for the in- 
applicability of the VT to interstellar clouds, since it will 
induce a flux of mass, moment and energy between the 
clouds and their environment, and will favor local collapse 
while disrupting the clouds within a dynamical timescale. 
The common assumptions discussed in the present contri- 
bution drive our understanding of the dynamical state of 
the interstellar clouds toward a picture that favors a static 
ISM. However, they are highly difiicufl to fulfill if the ISM is 
highly turbulent J as it was found to be many years ago (e.g., 
iMcCrav fc Snow^ll979^ . Inferences of the star formation ef- 
ficiency for supersonic (Mach numbers ~ 20 — 40) clouds in 
virial equilibrium living several dynamical timescales (e.g., 
Krumholz & McKcc 2005; Tan ot al. 2006) should be taken 
with caution. 

The lack of observational evidence for clouds being in 
VE, and the identificaion of asymmetrical line profiles ob- 
served toward interstellar clouds using different tracers are 
the best evidence of clouds being out-of-equilibrium sys- 
tems. These facts lead us to the conclusion that clouds 
should be transient structures, which exchange mass, mo- 
mentum and energy with their environment. This is precisely 
the opposite point of view to the old one in which clouds 
should be at rest, which still is present in textbooks and pa- 
pers, but it is consistent with recent observationally-based 
works that favor a scenario of rapid cloud and s tar forma- 
tion (see, e.g., the reviews bv iMac Low and Kle sscn 200^ 
iBallesteros-Paredes et al]l2006l and references therein.). 
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